IPPP/07/29 
DCPT/07/58 



Unparticles at heavy flavour scales: 
CP violating phenomena 

Roman Zwicky Q 

IPPP, Department of Physics, University of Durham, Durham DHl 3LE, UK 



Version of February 3, 2008 



Abstract: 

Coupling the scale invariant unparticle sector to flavour physics and assuming that it 
remains scale invariant we investigate its consequences in heavy flavour physics. A drastic 
feature of unparticle physics is an unusual phase leading to novel CP violating phenomena. 
We consider the CP asymmetry in the leptonic decay — > r"^z/ and the hadronic decay 
Bd — »■ D~^D~ , taking into account constraints of branching ratios and time dependent CP 
asymmetries. Generic plots are shown and it turns out that there exist parameters for 
which the CP violation is maximal. A prediction of a large CP asymmetry in B^ D^D^ 
is difficult to achieve in other models without contradicting the current data in other 
channels. The prediction of a CP asymmetry in leptonic decays, such as B^ r+z/, is 
novel. We identify the CP compensating mode due to the unparticles and show explicitly 
that it exactly cancels the CP asymmetry of B^ r^v as demanded by CPT invariance. 
Building up on earlier works we investigate the breaking of scale invaricance, due to 
the coupling to the Higgs and the size of the effects in the weak sector resorting to a 
dimensional analysis. An enhancement is observed on the grounds of the relevance of the 
unparticle interaction operator as compared to the weak four-Fermi term. 



*Roman. Zwicky® durham .ac.uk 



1 Introduction 



The possibility of a non-trivial scale invariant sector, weakly coupled to the Standard 
Model (SM), was advocated by Georgi in [1]. A scale invariant theory does not contain 
degrees of freedom with isolated masse^ unlike field theories used in phenomenological 
particle physics. Georgi called the degrees of freedom of such a theory " unparticles" . 

A non-trivial scale invariant sector, i.e. non-vanishing coupling, exhibits power-like 
scaling unlike the logarithmic scaling of QCD at the perturbative or trivial ultraviolet fixed 
point. The power-like scaling in Minkowski space seems to lead to curious phenomena. 
For example, the phase space of an unparticle with scaling dimension du-, which consists of 
the classical plus anomalous dimension, looks like a number of du (possibly non-integer) 
massless particles [T] . This could lead to interesting signals of missing energy. In a second 
paper [3] Georgi has pointed out that the unparticle propagator has an unusual phase 
^-iduir I sin((i2^7r) leading to spectacular interference patterns. 

By parametrizing a variety of interactions, unparticle phenomena were investigated at 
various energy scales and domains of particle physics such as electroweak physics [1], [S], 

collider physics [S], [S], [ID] (the latter investigates the (pseudo) resonance structure 
due to unparticles), DIS [II], [I2] 5, D-physics [I3], [H] [I5] [IE] [IT] , light flavour 
physics [18], [19] — 2 [13], [20] lepton flavour violation [21], [22], invisible decays [23] . 
cosmology [21] long-range interaction [25] and gravity [26]; All studies are based on the 
assumption that the theory remains scale invariant until the respective energy domain. 

Papers addressing questions of interpretation and the range of scale invariance have 
appeared. In an illuminating paper by Stephanov the continuous spectrum of the un- 
particle fields is discretized allowing for interpretation in terms of the language of particle 
physics. The authors of reference [5] address the question of the range of scale invariance. 
In the case where the unparticle couples to the Higgs vacuum expectation value (VEV) 
the latter will render the theory non scale invariant. This raises the question whether 
unparticle effects are observable in low energy experiments. A follow-up paper has ap- 
peared [6], where it is observed that higher dimensional operators can be parametrically 
enhanced under certain conditions on the scaling dimensions. Moreover this paper con- 
tains many physical applications and considers LEP-results to set bounds on the effective 
suppression scale. 

At first sight it seems rather difficult to pursue an analysis in low energy physics. The 
unparticle effects are parametrized in terms of an effective field theory where no principle 
is (yet) known to constrain the coefficients and a coupling to the Higgs VEV would take 
the theory away from scale invariance. On the other hand the novel phases could give rise 
to such striking phenomena that an investigation seems worthwhile. Moreover we have 
adapted the analysis of Ref. [5j to the weak sector and find that effects are possible if the 
unparticle field couples weakly to the Higgs VEV. The size of this coupling is not dictated 
by any principle and we may therefore regard its smallness as a working assumption. 

^In the absence of a mass scale the only possible particle candidates seem to be massless fields, but 
these have been shown to be free fields 0. 
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In the specific model or parametrization used, the unparticle will play the role of the 
VT-boson or charged Higgs in flavour-changing decays. It is well known that (time inde- 
pendent) CP violation manifests itself if there are at least two amplitudes with different 
relative strong (CP-even) and weak (CP-odd) phases. The phase of the unparticle prop- 
agator is CP-even and if we therefore allow for a different weak phase in the unparticle 
sector the door is opened to novel CP violating phenomena. Decays with one dominant 
weak amplitude seem particularly suitable, since they do not exhibit sizable CP violation. 
Moreover the unparticle should propagate at large heavy flavour energies because of the 
breaking of scale invariance. 

We analyze leptonic decays of the type B tu, where the SM and Beyond the Stan- 
dard Models (BSM) do not predict a CP asymmetry and Bd D~^D~ , which is further 
motivated by the unexpectedly large CP asymmetry measured by the Belle collabora- 
tion |43j. We shall also investigate how large the impact of unparticles can be without 
conflicting with branching ratio and indirect CP asymmetry predictions. 

The paper is organized as follows. In section [2] the scenario of the model is described 
including our parametrization of the effective Lagrangian and some general notation for 
CP violation is introduced. In section |3] the leptonic decay B ^ tu and B^ D^D^ 
are investigated followed by a discussion of similar channels. In section |4] we verify a 
constraint on CP-violation from CPT-invariance; namely that the partial sum of particle 
and antiparticle rates, with final states rescattering into each other, are equal. In section 
|5]we present the dimensional analysis of [5] adpated to a weak process. The paper ends 
with a summary and conclusions in section |6j 

In this paper we shall adopt Ay = 1 TeV as the scale of the IR fixed point. It is not 
difficult to rescale the results to a different scale, in the relevant places [Ay / 1 TeV) will 
be shown explicitly in the formulae. 

2 Scenario 

According to [1], with slighly adapted notations from [5], we shall imagine that at a 
very high energy scale ^ 1 TeV the particle world is described by the standard 
model fields and a self interacting ultraviolet sector. These two sectors interact with each 
other via heavy particles of mass Mu- The ultraviolet sector is supposed to contain a 
non-trivial infrared (IR) fixed point. An example mentioned in [1] is the Banks-Zaks [28] . 
perturbative type, fixed-point. Other examples are gauge theories with fermions in higher 
dimensional representations which exhibit near conformality, also known as walking. The 
phase diagram for an arbitrary number of fiavours and colours was given in [29] and 
preliminary lattice studies seem to confirm the theoretical expectations |30]. Below the 
scale My the theory may be described by non-renormalizable interactions 
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analogous to the four-Fermi interactions connecting the lepton and quark famihes within 
the SM. The ultraviolet theory flows into the IR flxed point around some scale Ku ~ 
1 TeV which will lead to new degrees of freedom called unparticles described by operators 
OiR = Ou- Below the scale Ku the theory may be described by an effective field theory 
in terms of the new degrees of freedom 

~ .d«+(lM-4)^SM0i, (2) 

with coupling A = cu{^u / MuY"^^^'''^'^''^'' and matching coefficient cu- Since we do not 
have a concrete theory at hand cu will be a free parameter to be constrained by experi- 
mental data. 

To make use of Lagrangians of the form ^ it will prove sufficient to know the cou- 
pling of an unparticle degree of freedom with given momentum P to its field operator 
for calculating decays into unparticles. Moreover via the optical theorem or dispersion 
relation it is possible to obtain the propagator and study interference effects, at the tree 
level, of unparticles and SM particles. In this paper we are concentrating on the latter 
effect. 

The propagator may be defined from its dispersion representation 
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^u{P'')^^j d'xe^^-{0\TOu{x)OlmO)= I -^^^^^^^ + s.t. (3) 
Jo JqT^s — r— IK) 

where we have assumed that the unparticle state \P) satisfies > and Pq > 0. 
The abbreviation s.t. stands for possible subtraction terms associated to a possible non 
convergence in the ultraviolet. The imaginary part is given by 

Im[A^(p2)] = |(0|O^(0)|P)p = A,,{Py--\ (4) 

The P^ dependence is solely determined by the scaling property of the operator Oy- The 
factor Adi^ is a priori an arbitrary normalization constant which has been chosen to be 

in reference pp. It is the phase space volume of du massless particles. This choice was 
motivated by the fact that Q exhibits the same functional behaviour as du massless 
particles. This together with the fact that unparticles would (presumably) escape particle 
detectors has led Georgi to point out [1]: "Unparticle stuff with scale dimension du looks 
like a non-integral number du of invisible particles." With Q the dispersive integral ^ 
is elementary: 

' 2sin(di^7r) (-P2 -i0)2-rf« 2^m{duT^) [P'^f-^u ' ^> 

We observe that the unusual phase, due to analtytic continuation to Minkowski space, is 
due to the non- integral scaling dimension du- It has been shown in [3] that the discon- 
tinuity of the propagator yields the imaginary part (|4]), which we have already implic- 
itly assumed in the dispersion representation ([3]). The phase is accompanied by a term 



ds Im[Aw(s)] 
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sinlduir) in the denominator. In simple cases of direct CP asymmetries this term will 
cancel and in more complicated cases it will partly cancel. The cancelation of this term 
will play a crucial role when we verify the the equality of partial rates of particles and 
antiparticles, which is a consequence of CPT, in section [4| Plots of the function A^^ and 
Adi^/ sm{di(7i) can be found in the appendix [A} 

A strong (CP-even) phase in a propagator appears spectacular and is not an element 
of common models in particle physics. The phase factor is due to dynamics in the un- 
particle or scale invariant sector. The exactly solvable two-dimensional Thirring model, 
which contains fermions with a current-current interaction term, is an example where the 



dynamical phase can be seen explicitly, c.f. appendix B.2 The anomalous dimension is 
a function of the coupling constant and assumes the free field value in the case where 
the coupling is taken to zero. We would like to add that in the Thirring model, due 
to the fermion selection rule, the anomalous phase is not so immediately observable. It 
is also interesting to note that there is a connection between the non-trivial phase and 
the causality structure. The commutator of the unparticle field in the vacuum is, c.f. 
appendix B.l Eq. ( A.6[ ), 



(O|[0^(a;), 0^,(0)] |0) = -t sign{xo)e {x') {x')-'"'^ /(du) , (7) 



where f{dn) is a function explicitly given in (A. 6), The commutator vanishes for space-like 



X < separation and obeys causality. For generic du the support is inside the light cone 
which seems in agreement with the spectrum > 0. For an integer value, dy = n, the 
commutator behaves as ~ S^"~^\x'^) and has support on the light cone only. Note that 
in the latter case there is no CP-even phase in the unparticle propagator. 

We shall now discuss the possible values of du- In the upper range du = 2 is singled 
out, since we observe that for du < 2 no subtraction terms are needed. Moreover the 
singularity of du approaching an integer value larger than 2 has been interpreted in p] as 
the a du-psxticle cut which should not be attempted to be described by a single unparticle 
field. In the lower range the value du = 1 seems special since it corresponds to the free 
massless field, 

limA^(0 = 1. (8) 

It is also observed that for du I the dispersion integral does not converge in the infrared. 
This might be interpreted by the fact that the field decreases even slower than the free 
massless field in coordinate space. It has been shown that for dy < 1 the conformal group 
does not admit unitary representations [32]. Moreover in reference |1] it was noted that 
the decay into an unparticle has a non-integrable singularity in the decay rate for du < 1- 
We shall therefore think of du as being 

l<du<2 (9) 

or parametrising with repsect to the free field limit, du = l + 'j with anamalous dimension 
< 7 < 1 at the non-trivial fixed point. 
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2.1 Parametrization of the effective Lagrangian 

In this section we shall give our parametrisation of the coupling of the unparticle sector 
to the SM. We will investigate charged-flavour decays and therefore it is sufficient to give 
the couplings to that sector. We couple a vectorial unparticle operator Ou to a scalar 
and a pseudoscalar density. The unparticle will therefore be a charged Lorentz-scalar and 
play the role of a charged Higgs rather than a PF-boson. We parametrize the effective 
Lagrangian as follow^ 

+ #i(?"75g) Ou + -J^(^750 Ou + h.c. , (10) 

where q' = {u,c,t), q = {d,s,b), v = (i^e, t'/^, i^t) and / = {e,fi,T) are summations over 
the families. In the notation of Eq. (|2]), dsu = 3. The weak (CP-odd) phases are 
parametrized as deviation from the phases of the CKM matrix V^/q and analogously the 
leptons as deviations from the PMNS matrix 11^1. 

= e'<'.\Xi^\ = arg[l^,g + <50j^ . (11) 

The Lagrangian is a non vectorial copy of the charged current sector in the SM. This 
allows us to apply up to some level the same tools in the unparticle sector as in the SM. 
Note that the unparticle carries charge, unlike the Lagrangian used in [Tj. 



We would like to stress that the Lagrangian in Eq.(lO), in the absence of an explicit 
model realising the unparticle scenario, is not dictated by any structure and is therefore 
only an example. Other Dirac structures with the same flavour transitions are possible. 
The axial and vector structures, for example, can be coupled to a transversal unparticle 
Ojj i^f^Ou = 0) or to a derivative coupling of a scalar unparticle d^Ou 

SC^^ = ^S^(9"7.(75)g) d^Ou + ... (12) 

The former leads to a propagator transversal propagator ~ {—Qfiv + PiiPv/p^) which van- 
ishes by transversality when coupled to a pseudoscalar particle of momentum p^, which 
is the case in the examples considered in this paper. The latter leads to an identi- 
cal contribution as the scalar and pseudoscalar contribution in the examples considered, 
^s,p ^ const ■ \d{v,A), where Xq(v,A) is parametrically supressed by one power of Au as 
compared to Xs^p- We shall give the results in terms of Xs^p in the paper but also indi- 
cate explicitly how they change for a A9(y_A)-coupling. Please note that although we have 



^The channel i?+ r+i^ is mediated by a {PxS) + (PxP) structure whereas D~^D^ decays 

via a {SxP) interaction. The vector and axial couplings are discussed in the text above and we do not 
consider tensor couplings since they do not couple to single scalar particles. 
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just stated that Xs,p and the Xd[v,A) equivalent in the examples considered, their role 
in model building, in regard to SU{2)l for example, might be rather different since the 
former couples fermions of opposite chirality whereas the latter couples fermions of the 
same chirality. 



2.2 General formulae for branching ratios and CP asymmetries 

In this section we shall give the formulae for the branching ratio and CP asymmetries, 
used later on, in the case of two amplitudes with different strong and weak phases. This 
paragraph is completely general in principle, but we shall have in mind that one amplitude 
is due the SM and the other is due to the unparticles. For a decay B CD we parametrize 

A{B ^X) = Aie'^^e'^^ + Aic'^^e'^^ , (13) 

where the 5j denote the CP-even phases and (pi the CP-odd phases. The branching ratio 
B and the CP averaged branching ratios are given by 



where 



B = B'U B = B'^U, (14) 



B' = T{B)^^^Xy'{ml,ml,ml)\A,\' A 



/a = (l + 2Acos(</.i2 + 5i2) + A') 

/a = (1 + 2A cos(0i2) cos(5i2) + A^) , (15) 

and A(a, 6, c) = + 6^ + — 2{ab + ac + be), 0i2 = 0i — 02 and S12 = 5i — ^2- In the case 
where the transitions 

B ^CD ^ B (16) 

are possible, the 5-meson is neutral and the dynamical mixing of B^ and B^ leads to 
a time dependence in the CP asymmetry. In the case where the coefficients q and p, 
relating the flavour and mass eigenstates of the neutral system, assume \q/p\ = 1, the 
lifetime difference AT/ AM ^ 1 the CP asymmetry assumes the following form 

AcAB, - CD) ^ If^ZcDj+vlBZcD] = M^Mt) - Ceo cos(AMt) . (17) 

Both assumptions mentioned above are satisfied for the Bd system. The sign convention 
is such that the b ^ cd decay rate enters with a plus sign, please note B = Bq = Bd 
ibd) [35J. Writing q/p = e~^'^'^, where ipd = 2(3 is the mixing phase of the Bd system 
A = q/p{A/A), then the (time independent) CP asymmetry assumes the following form 

1 _ I \|2 2A 

C=—\T\-,=^sm{5,2)smi^,2) • (18; 
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In the case where the system CD is a CP eigenstate with eigenvalue ^cd = ±1 , which is 
a particular realization of (16), the time dependent CP asymmetry assumes the following 
form 

^ = ^CD f ^"'i^il = ecD9^(sin(0rf + 20i) + 2Acos(5i2) sin(0d + 0i2) + sm{(t)d + 2^2)) • 

J- + I'^l /A 

(19) 

For Bd J/'^Kg ij/<i!Ks = —1, 01 — and there is no sizable second amplitude in the 
SM and therefore A ~ and 02 and the gold plated formula Sj/^ik^ = sin(2/3) is 
recovered. 



3 A leptonic and a hadronic decay 
3.1 r'^u; scale invariant sector at 5GeV 

In the standard model charged pseudoscalars decaying to a lepton and a neutrino are of 
particular interest because of their simple dependence on the pseudoscalar decay constant 
and the CKM matrix element, see below. 

The novel feature when adding unparticles is a CP asymmetry. We will investigate 
how large this asymmetry can be, remaining consistent with the the branching ratio 
measurement 

Below we will give the decay amplitude for r^u for the SM and the unparticle 



contribution with effective Lagrangian as given in (10). 

The unparticle is propagating at the scale and we therefore assume that the scale 
invariant sector extends down to the scale. The SM and unparticle graphs are shown 
in Fig. [T| where we have indicated the phase. 

The additional unparticle amplitude leads to a slight complication. As a matter of 
fact in experiment we do not observe the neutrino flavour but an inclusive measurement 
on the neutrino flavour is performed since the neutrinos are not detectable. In the case 
where there is only one amplitude, as in the SM, unitarity of the PMNS matrix hides 
this fact from the final formula. This is not the case for the unparticle amplitude and we 
shall therefore derive formulae for B~^ t^u = 'Y^iB^ r^i^i via B^ r^i^i- The 
amplitude is the sum of two incoherent terms of opposite parity in the final stat^ 

A{B+ ^ = ^KTb^r., fsmr ■ 

([Pr](l + ^l^e-'^^-e-'^') + \vi,t]{-1 + Af^^e'^'^^-e-^^r )) , (20) 



''The amplitude (20 1 displays the famous hclicity supression in the SM due to its chiral structure which 
manifests itself in the fact that the amplitude is proportional to the lepton mass. For a pseudoscalar 
coupling, as in the charged Higgs model, or the one used here, the helicity supression is relieved as can 



be inferred from Eq. (21 ). If we were to use a derivative coupling to the axial vector \12\ instead of 
a pseudoscalar coupling as in (10) then the following substitution, ui^/{mi,mr) rn^/Ali, in Eq. (21 1 



would reproduce the result for the derivative coupling. 
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where = Scp^^ — Scp!^^^ for D = {S, P), I = (e, /x, r). The i?-meson decay constant is 
defined as rrih {0\bi'y5u\B~^) = fB^\ , where we neglect isospin breaking effects. The ratio 
of unparticle to SM amphtude is 

I \'^'^' I 

^TUl — ITT ,^TU—ri L^TV 1 

\'Jtvi I 

\\f\ A,^ ml {Gf/V2)-^ (mlY^-^ 



|Vub| 2sin{duT^) rrii^mr ^ ^ 

We will now make a simplifying assumption in order to simplify the analyis. We impose 
the left-handed chirality on the unparticle sector i.e. A^*^' = —X^P and {^rv = Ai^iJ^\ 
6(j)ru = fi = rp'^"*). This means that the amplitudes for opposite parity give the 

same result and this allows us to combine the two amplitudes into one. The branching 
fractions to a specific neutrino fiavour final state are 

B{B^ - r^u,) = Blf Va., B{B^ - r^u,) = B^f^ |t/..j7A., , (22) 



with / and / as in (15), 0i2 = — 0z, 5i2 = —du7r. The familiar SM branching fraction 
reads 

= r{B^) # K.\'fBmsml{l - "^f (23) 



and does not depend on the neutrino fiavour. Please note that in the SM B^^ = B^ff^. 
The experimentally tractable or neutrino inclusive branching fraction is 



B{B+ r+z/) = J2B{B+ r+ui) 
I 

= Blf J2 \Ur.^\l + 2nA., cos(0, + dun) + (nA.,)^) 
I 

= Bf^il + \Uru,\\2riAr,cos{<l)i + du^) + (riA^u?) ■ (24) 

The formula could be further simplified if the ri were independent of /, which we shall 
assume shortly below. The CP averaged branching fraction is 

B{B^ ^ T+u) = B^^J^ = B^,^il + J2\Uru^\2rlK.cosi(Pl)cosidu7r) + irlAr^^^^^ (25) 
The CP asymmetry assumes the following form 

, , , T(B- ^ T-U) -T(B+ ^ T+U) 2Aru . ^ , ,9 . ^ 



where JF is implicitly defined in Eq. (25). Let us note that the CP violation encoun- 
tered here is proportional to ~ lm[V*^\^Uru^s^*], which is hidden in the formula above, 
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Figure 1: (left) SM diagram for B tu (riglit) unparticle diagram with CP odd phase e^'^^'" . 
The unparticle is denoted by a double line. 

and is the product of two quadratic reparametrization invariants. The effect is entirely 
proportional to the sine of the phase difference between the CKM (PMNS) and the un- 
particle flavour sector and can therefore not occur in the SM. In order to do a qualitative 
assessment we shall study the case where there is no flavour dependent perturbation in 
the neutrino sector and therefore drop the label /. The formulae for the CP averaged 
branching ratio and the CP asymmetry then simplify to 



T V] 



Blfil + cos((/)) cos{duTi) + A 



=±7r/2 



TV / 



ll^rv sin(0) sin((ii^7r) 



2A^^ cos(0) cos{duTT) + A2^ 



=±7r/2 ±2|A^i,|| sin((ii^7r) 
^ 1 + A2 



(27) 



where in the last step we have simplified the formulae further by setting the weak phase 
difference to 90(270)°|5 N.B. in the notation used in Eq. Acp{Tiy) = -Cry. This 
choice maximizes the CP violation for appropriate values for Aj,;. Before we are able to 
constrain the CP violation with the rate we have to give the theoretical and experimental 
results of the latter. 

The following hadronic parameters, t^^ = 1.643 ps, fs = (189 ± 27) MeV a lattice 
average from [35] and \ Vub\ = 3.64(24) ■ 10~^ from the fit to the angles of the CKM triangle 
[36], are used to estimate the SM branching fraction 



T V 



SM 

theory 



83(40) ■ 10" 



(2^ 



We have doubled the uncertainty due to |Vub|- This estimate has to be compared with 
the measurements at the 5-factories 



units 10"^ 


i3(5+ ^ r+z/) 


BaBar[37](223MBB) 


90(60)(10) 


Belle[38](449MBB) 


179(53) (48) 


HFAG[3l] 


132(49) . 



(29) 



^N.B. sin((i;^7r) < for 1 < rf;^ < 2 as assumed throughout this paper (|9|. This is the reason for the 
absolute values in the equation above. 
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Weak phase = 90(270)°, flavour independent perturbation neutrino sector 



B(B*^ r*v) [10"''] 




C (B*^ i:*v) [ I sin (du/r) 
1 




Figure 2: A weak phase difference (/> = 90(270)° is assumed here for A,y^ positive(negative). 
(left) Branching fraction (27) as a function of The black bands correspond to the SM 



estimate (28) at A^,y = 0. The blue band corresponds to the HFAG bounds in Eq. ( |29| ). (right) 
The CP asymmetry as a function of in units of | sin(dz^7r)|. The scale k.u = 1 TeV is chosen 
here. N.B. in the notation used in Eq. (26) Aqv{j'>^) = —Cry 



In Fig. [2] (left) the branching fraction (27) is plotted as a function of l^^v with un- 

= 0. The blue band corresponds to 



certainty taken from the SM estimate (28) at A 



the HFAG bounds in Eq. (29). The CP asymmetry is plotted to the right of that figure. 



The branching ratio does not set limits on the amount of CP violation, demanding the 
uncertainty bands to be tangent at worst |A^;^| < 1.8. Even in the case where the HFAG 
and theory uncertainty are halved, the value |A^i,| = 1, at which the CP asymmetry is 
maximal, is still consistent. 



Weak phase </> ^ 90(270)°, flavour independent perturbation neutrino sector 

In this subsection we shall repeat the analysis for a general weak phase difference and 
show two dimensional plots in the variables (0, du) for different ratios of effective couplings. 



The quantity A^-j, (21), used in the previous paragraph, depends on the ratio of effective 



coupling and scaling dimension as follows 

A^^ ml {Gf/V2)-' 



A 



Tl/ 



m 



B 



where 



^'^ 2 sin{dij7T) 1711,171^- ml ^ ^ti 

Pr. ^ ^TTX • (30) 

A plot relating A^j^ and dy can be found in appendix |A| Fig. |9} In Fig. [s] (right) the 
the CP asymmetry Cru is plotted as a function of {(p,dii) for p^-^, = (10°, 10~^, 10~^). 
The pattern is clearly regular and the condition for a large asymmetry is \Aru\ ~ 1- 
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For smaller values of pru the amount of possible CP violation is decreasing because the 
condition mentioned above cannot be satisfied. The constraint on the branching fraction, 
Fig. |3] (left), is defined by the following acceptance function 

A{du, 0, p) = (1 - r{du, 0, p))6(l - r{du, 0, p)) , 

r{duA,p) = ^|C(l + 2A,,cos(0)cos(rf^7r) + A^J-S^^^G| ^g^^ 



for = 83 ■ 10-^ ijHFAG ^ . lo-^and for the quantity A5 we add the uncertainty of 
the SM prediction and the HFAG value linearly to AS ~ 80 ■ 10~^. This function assumes 
values between and 1, where 1 signifies maximal agreement and means that the point 
is excluded; or in other words we consider predictions with a deviation larger than /S.B 
as excluded. 



For smaller values of the linear term for the branching ratio in Eq. (27) becomes 
dominant and a regular pattern in cos(0) emerges. Note that since the predicted branching 
fraction is lower than the central value from experiment, the weak angle = 180° is 
currently disfavoured since it would lower the theory prediction even more. 

3.2 Discussion and remarks on B fiu, D fiu, Bg fi'^fi^ etc 

We have seen that applying the unparticle scenario to the leptonic decay B ^ tv leads 
to CP violation. There is no experimental data available that gives both the negative and 
positive charged semileptonic decay rates, i.e. quotes (bounds) on CP asymmetry in a 
semileptonic decay. 

The current data on B ^ tv do not allow us to set bounds on the amount of pos- 
sible CP violation. The amount of events at BaBar and Belle are of the order ~ 20. 
An improvement in theory, in particular on the 5-meson decay constant, and the large 
statistics of a Super 5-factory would of course improve the situation. Unfortunately the 
decay B^ ^ t^v will not be possible or competitve at LHCb because of the neutrino final 
state and the intricacies in the r detection, whether D{Ds) {t, p)^ decays are possible 
at LHCb is currently under investigation. 

We shall comment on other leptonic modes. They are all described by the same 



formula (23) for B ru with obvious substitutions for V^b, /s, ms and m^. We may 
also consider the D-decays assuming that the scale invariant sector extends to ~ 2 GeV. 
The decay p'^v is measured by CLEO [32] the ~ 50 events lead to a thirteen 

percent accuracy. The decay constant = 220(20) MeV is taken as an average value of 
theory determinations from the table in [39] and \ Vcd\ = 0.227 [35]. The Cabibbo allowed 
decays p^u are measured as well [35], although with less precission. The decay 

constant /^^^ = 264(36) is obtained from an average of f^Jf^ = 1.20(5) of the table 
in [22] and \Vcs\ = 0.957(17) (93) [33]. A summary of the experimental [53] and theory 
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Figure 3: A horizontal line of figures corresponds to different fractions of effective couplings 
p = lo'^'^^'^^ as defined in (30). (left) Constraints on the {(p^du) parameter-space from the 



branching fraction. The values in the dark regions are allowed whereas white ones are excluded, 



c.f. Eq. 31 for a more details, (right) The CP asymmetry as a function of {(jy^du)- The scale 
All = 1 TeV is chosen here. 
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predictions is: 





B ^TU 


B ^ fiu 


B ^eu 


Experiment 
Theory 


132(49) • 10^6 
83(50%) ■ 10-6 


< 17-10-^ 
3.7(50%) ■ 10-^ 


< 9.8 ■ 10-*^ 
8.4(50%) ■ 10-12 




D ^TU 


D ^ fiu 


D ^eu 


Experiment 
Theory 


< 2.1 • 10-3 
1.1(20%) ■ 10-3 


4.4(7) • 10-4 
4.3(20%) ■ 10-4 


< 2.4 ■ 10-5 
1.0(20%) ■ 10-10 






Ds fiu 


Ds — * eu 


Experiment 
Theory 


6.4(15) • 10-2 
5.5(30%) • 10-2 


6.3(18) ■ 10-3 
5.7(30%) ■ 10-3 


not available 
1.3(30%) ■ 10-^ 



The i? decays are predicted to 50% due to uncertainties in and |Vub|, whereas the 
D{Ds) decays have a lower uncertainty 20(30)% due to foifDs)- The helicity supression 
in the SM is apparent from the table. 

Repeating the analysis for D~^ Z^^^^, as shown in Fig. |2| we obtain that \Ad^hi^\ < 
0.65 which still allows for a rather large CP asymmerty, \Cd^iiu\ < 0.9. 

The prediction of these modes in the SM is solid and a significant deviation would 
be a clear hint for new physics. In particular one expects larger rates in models where 
the helicity suppression is relieved. An example is the charged Higgs or the effective 
Lagrangian used in this paper. The charged Higgs does not predict a significant CP 
asymmetry whereas in unparticle models it is possible and therefore a CP asymmetry 
could be used to discriminate between the models. 

We would also like to mention the decay /^^^, the KLOE collaboration reports 

~ 860 events and a branching ratio B{K^ /i+i/(7)) = 0.6366(9)(15) |10]. On the one 
hand it seems unreasonable that the scale invariant sector could extend to ~ 500 MeV 
but on the other this channel has the largest statistics. If we assume that theory predicts 
the rate to 5%(10%) this would roughly bound \/S.k^^u\ < 20(30) and the CP asymmetry 
to \Ck^^,\ < 0.4(0.55). 

Finally a comment about B(d,s) This channel is rare since it is a flavour- 

changing neutral decay further suppressed by the coupling of the Z and the helicity of 
final states, E{B(^cl,s) fi~^fJ'~)^^ ~ 10-iO(10-*). The branching ratio is not yet measured, 
the bounds are about one and half order of magnitude away from the SM prediction. An 
analysis along the lines of B —>■ ru does not make sense since there are no direct constraints 
in that channel. A possibility would be to combine it with constraints from AM^a^s), which 
are measured, as advocated in reference [H]. 
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3.3 Bd D^D ; scale invariant sector at 2 GeV 

The decay — D^D^ corresponds to a 6 — ccd transition at the quark level and is 
colour allowed. It has the same quark level transition as J/'^ttq but two complica- 

tions arise as compared to the latter. First, since it is colour allowed it receives sizable 
contributions from a gluonic penguin [H] and second the final states combine into a sum 
of isospin / = and 7 = 1 waves which have in general different final state interaction 
phases. Ultimately we will neglect the penguins in our analysis, to be discussed below. 
Our motivation to investigate the Bd D^D~ is driven by the measurement of a large 
CP asymmerty by the Belle collaboration [43 The SM expectation is C|,'f^_ ~ —0.05. 





Cd+d- 


Sd+d- 


BaBar[42j(364MBB) 


0.11(22)(07) 


-0.54(34)(06) 


Belle |13](535MBB) 


-0.91(23)(06) 


-1.13(37)(09) 


HFAG 


-0.37(17) 


-0.75(26) 



(32) 



It has to be said that the Belle result is somewhat moderated by a significantly lower 
value from BaBar [12] with opposite sign. Note that the cental values from Belle also 
violate the general bound + 5^ < 1. 

It shall be our goal to see how large a CP asymmetry Cd+d- the unparticles scenario 
can generate and still be consistent with the branching fraction and the time dependent 
CP asymmetry. 

In our analysis the unparticle will replace the W in the tree level amplitude in, c.f. 
Fig. |4] (left). We therefore assume that the scale invariant sector extends to the D-meson 
scale ~ 2 GeV. 




Figure 4: 5 ^ dec (left) tree diagram, (right) penguin diagram 

We shall first reconsider the situation in the SM before we move on to the unparticles. 
Writing the amplitude as the sum of the tree and penguin topology 

A{Bd D+D-) = At + Ap = At{1 - e'^^^e'^r^^) , (33) 
am grateful to Christopher Smith for drawing my attention to this measurement. 
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the ratio of penguin to tree amplitude rpT can then be estimated by the Bander-Silverman- 
Soni mechanism [JT], c.f. [H] or |l5] for an updated analysis, 



^PT 



0.08 5pT^ 205° . (34) 



This allows us to obtain the asymmetries from (18) and (19) 



C|V^-0.05 5|V^-0.78 . (35) 



Comparing with the experimental results (32) we infer that the SM is in good agreement 



with the time dependent CP asymmetry Sd+d-- The direct CP asymmetry Cd+d- — 0.05 
is about two standard deviations lower than the HFAG value 0.37(17). In view of the 
non consistency of the two measurements it is certainly wise to wait for updates from 
the i?-factories. We will in the following neglect the penguin contribution in regard to 



its moderate size (34) in the SM. We will also neglect the "unparticle penguin". The 
ratio of the unparticle penguin amplitude to the unparticle amphtude is expected to be 
of the same size as in the SM, unless the up-type transition is enhanced by the effective 
couplings. We are therefore implicitly assuming that l^^^^p^A^^p^l ^ I -^(5 p) -^(5 p) I • 

We will describe the amplitude Bd D^D^ within the naive factorization approxima- 
tion. Naive factorization describes colour allowed modes (topology as in Fig. |4]to the left) 
like B 7r+7r+ and Bd D^n^ with at least one fast or light meson with an accuracy of 
around 10 — 20% level. For Bd D^D^ , factorization in general and naive factorization 
are not expected to hold. The overlap of the emitted D+-meson with the Bd 
transition is expected to be relatively large. However it is empirically observed that naive 
factorization still works reasonably well. We shall account for final state interactions, not 
included in naive factorisation, by an isospin analysis which is presented in the appendix 
[Cj The effect is that the amplitude receives a contribution cos((^i — Sq)/2) ~ ±0.63(15), 



c.f. (A. 17). In fact the sign is not determined but since it enters in the square in the 



observables it is of no concern here. The amplitude for Bd — > D^D in the SM is 

A{Bd - D^D-) = ^K;Kdai/D((m| - ml)f^^{ml) + m^/^^(m^)) 

X cos((5i-5o)/2)e^(^-^«)/2 = ^SM^ 

where ai = C2+C1/3 ~ 1 is the colour allowed combination of tree level Wilson coefficients 
and the D-meson decay constant is defined as nic {0\ci'y5d\D~) = fofni), where we neglect 
effects due to isospin breaking. The B ^ D form factor can be parametrized by use of 
Lorentz covariance as 

{D\h,c\B) = f^''{q'){pB+PD), + f^''{q')q,, (37) 

with momentum transfer q = Pb — Pd- The form factors are related to the famous Isgur- 
Wise function /f^(g2) = ^^^^(^), /BD(^2) ^ -y/^^e(^) m the heavy quark 
limit. Here w = v ■ v' = {m\ + m\) — q"^)/ {217131710) ■ Whereas the normalization of 
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the Isgur-Wise function ^(1) = 1 follows from charge normalziation in the heavy quark 
limit the values around maximum recoil are much less known. We shall take the value 

= 0-54 from |17] and scale it up to = m|, by use of a single pole model 
^{w) ~ y^2/(w + l)(wmax— w(m^*)/(w— t(7(m^,). The i?*-meson has the correct quantum 
numbers = and its mass is the same in the heavy quark limit as = 6.29 GeV 



We obtain 



[m 



D) 



0.7. With fo = 220 MeV. we get 
B{B, ^ D+D-)l^,,,^ = 1.7(10) • 10-^ (38) 
as a theory estimate, where the bulk of the uncertainty quoted is due to the isospin final 



state interaction phases (A.17). This estimate has to be compared to the experimental 
value [55] 

i3(5, D+D-)pDG = 1.9(6) ■ 10-^ . (39) 



The agreement seems accidentally good in regard to the approximations made. 
As in the previous section we parametrize the amplitude 



(40) 



with J\^^ as given in 
unparticle amplitude i: 



and relative weak phase (f) = 5(j)ch — 5(j)cd- The ratio of SM to 



A 



DD 



1^5 I 1 



Ad,. 



m 



D 



m 



D 



m 



D 



A2 



(41) 



Note that, unlike for B tu, the negative parity of the D-meson selects only the A 



cd 



coupling in the final vertex. The observables are obtained from Eq. (18) and (19) with 
Cd+d- = 1, 0d = 2/5, 01 = 0, 02 = -0 and 612 = du-n- : 



B 



DD 



KjSM , 
■ UdD i Add 



B 



DD 



—(j) and 612 



SM 

DD i Add 5 



2A 



DD 



DD 



DD 



/Add 

-1 
/Add 



sin[0] sin[(i2^7r] , 
(sin[2/3] + 2 Add cos[du7r] sin[2/3 - 0] + A^^ sin[2/3 - 2(f)]) 



and 



^DD 



r(5. 



Weak phase 



327rms 
= 90(270)° 



im<-rnl)f:^iml)+mlf^^{ml)Y\V:M 



(42) 



(43) 



In order to look for maximal CP violation we may again set the weak phase difference 
to 90(270)° in the formulae in Eq. (42). In Fig. [s] (left) the branching fraction is plotted 



derivative coupling d^Ou to the vector Dirac-structure (|T2|) instead of a scalar coupling as in (flO' 



would lead to a change of mjj/{mc{mb — rric)) mj^/A^ in Eq. (41 1. 
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A 



Figure 5: A weak phase difference (f) = 90(270)° is assumed liere for Add positive (negative), 
(left) Branching fraction (42) as a function of Add- The black bands correspond to the SM 



estimate (38) at Add 



0. The brown- red band corresponds to the HFAG bounds in Eq. (39). 

as a function of Add for du 



(middle) Time dependent CP asymmetry Sj^+d- as a function of Add for dy = 1.1,1.5,1.S; 
where the dashes get shorter for larger values of du. The interpolation between those values is 
fairly smooth, (right) The CP asymmetry as a function of Add in units of | sm(di(TT)\. 



as a function of Add with uncertainty taken from the SM estimate (38) at Add = 0. 
The brown-red band corresponds to the HFAG bounds in Eq. (39). The new feature as 
compared to the B ^ tu analysis is the constraint from Su+d- which corresponds to 
the figure in the middle. The CP asymmetry is plotted to the right of that figure. Once 
more the branching ratio does not set limits on the amounts of CP violation, in fact the 
uncertainties are very similar as in i? ^ rv. Demanding the uncertainty bands to be 
tangential at worst results in |Add| < 1-5. The constraints from Sd+d- do depend on 
the scaling dimension. The parameter du = 1-1 for example seems slightly disfavoured as 
compared to the value djj = 1.9 



Weak phase ^ 90(270)° 

We investigate the two dimensional parameter space (0, d^) for different ratios of effective 



couplings. These quantities relate to Add (41) as follows 



^DD — Pdd; 



' 2 sin{du7!-) mdmb - m^) ml V A^ 

17.10»(3.5.10-«-4J^)*-'^Poo, 

V ITeV/ sm(aiY7r) 

I \ cb \ cd I 
I »/cb(^cd| 

A plot relating Add and du can be found in appendix lAl Fig. [91 In Fig. |6] (right) CP 



where 



asymmetry Cd+d- is plotted as a function {(t>,du) for pdd = (10°"^"''). The pattern 
is very similar in its form to B —>■ tv. A large asymmetry is obtained for |A^j,| ~ 
1, which cannot be attained for smaller values Pdd- The constraint on the branching 
fraction, Fig. |3] (left), and the CP asymmetry S £>+£,- are evaluated with the same kind of 



acceptance function as for i? rz/ (31 ). The corresponding values for the CP asymmetry 
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C(B^ D*D") 0=10" 




Figure 6: The observables with fractions of effective couplings p = 10 ' ' , as defined in (44) 



are plotted from the top of the figure to the bottom. Constraints on the ((/>, du) parameter-space 
from (left) the branching fraction and (middle) the CP asymmetry S £,+£>- (middle). The values 
in the dark regions are allowed whereas white ones are disfavoured, c .f. text for more details, 
(right) The CP asymmetry C^+d- as a function of {(p,di{). The scale Ak = 1 TeV is chosen 
here. 
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are S^j^^ = - sm(2/5) = 0.69, S^^^^ = -0.75 and = 0.52 corresponds to two 
standard deviations. The values for the branching fraction are i3|J^^_ = 1.7 ■ lO"'^, 
B^^. = 1.9 • 10~^ and AB = 1.6 corresponds to hnear addition of the theoretical and 
experimental uncertainty. 

A qualitative result that can be inferred from Fig. [3] is that the parameter space of a 
large positive CP asymmetry Cn+c- is disfavoured by the bounds from the Sd+d-- This 



is easily seen from the formulae (42), (41) and the plots in appendix (|Aj). A negative 
C £,+£,- demands a weak phase (f) < 180° and then the linear and quadratic terms in 
Sd+d-^ add constructively and are in conflict with the consistent result between the SM 
and experiment in this observable. As for B ^ tv for small pdd the linear terms dominate 
the quadratic ones and a regular pattern in cos(0) and sin(2/? — 0) emerges. 

3.4 Discussion of D^D and remarks on U-spin &; colour 

related channels 

A large CP asymmetry C^+d- would be a rather puzzling fact, as for instance discussed 
in Ref . [IE] . One is lead to suspect that the gluonic penguin Bd Dqq with q = c might 
be enhanced by new physics. This scenario would or should lead to enhanced penguin 
amplitudes for q = [u, d, s) as well and enter B^ (vttt, KK) in disagreement with the 
i?-factory data. 

We have seen that an unparticle scenario can lead, for appropriate parameters, to 
enhanced CP violation. One might wonder whether similar results shouldn't also show 
up in U-spin (s ^ d) and colour related channels. The plots in Fig. |6] indicate that the 
CP asymmertry S in general does not necessarily receive large contributions. This can 



be inferred from Eq. (42) or by noting that the unparticles just contribute to a large 
SM background from sin(2/3). We shall therefore focus on the CP asymmerty C. Let us 
note however that the situation for Bs decays is different since the mixing phase (ps — 
~ 2(3) in the SM and the contributions of unparticles would be not be shielded by a 
large SM value. 

The colour related or colour suppressed channel of B^ —>■ D^D^ is B —>■ J/^ttq. The 
CP asymmetry has been measured Cj/^t^q = —0.11(20) [35j, which is not conclusive in 
regard to its size. In the colour suppressed modes the non-factorizable contributions are 
enhanced due to different combinations of Wilson coefficients (typically ~ 2 — 3 larger 
than the factorizable amplitude) and have large strong phases. On the practical side it 
is harder to estimate them reliably in the SM and even more in the unparticle scenario, 
where the unparticle is dynamical as compared to the contracted VT-boson propagator 
in the SM. The strong phases and the different hierarchy between factorizable and non- 
factorizable contributions in the SM and the unparticle scenaricj^ make it impossible to 
draw conclusions without explicit calculations. 



parametric estimate gives that the non-factorizable contributions in the unparticle scenario are 
suppressed by a factor 2m|j/(mj^^ + m^) ~ 0.2 as compared to the SM. 



19 



The U-spin related transitions b ccs are CKM enhanced and therefore statistics 
should make them more attractive. In principle there is no reason that generic new physics 
respects the CKM hierarchy and U-spin. In the unparticle scenario there is no principle 
that dictates a CKM-like hierarchy in the coefficients Xg'q in the effective Lagrangian 
(10). Therefore they are not necessarily of major concern. Let us nevertheless discuss 
them. The gold plated decay B J/'^Ks is also colour-suppressed. The measurement 
of the CP asymmetry Sj/^Ks = sin(2/5) has allowed determination of the angle /3 in the 
SM, whereas the CP asymmetry Cj/xfiK^ = is consistent with experiment. This mode 
is highly consistent with the SM or more precisely with one dominant amplitude. The 
branching fraction of the colour allowed decay Bd D^D~ has been measured but no 
CP asymmetry has been reported, presumably because it does not exhibit CP violation in 
mixing. If the Belle CP asymmetry in Co+d- gets conffimed a look at the CP asymmetry 
appears mandatory. 

In summary the most interesting parallel channel is probably Bd J/'^ttO and the 
improvement of the measurement in Cj/^ttq should be watched along with C £,+£>-. In the 
scenario we described we would generically expect a large CP asymmetry Cd+d- to be 
accompanied by a large asymmetry in Cj/.^7ro- It is a serious point of criticism, but on 
the other hand the experimental result is not conclusive and in theory there might be 
cancellations between the strong phase e*'^"'^ and the phase from the non-factorizable in- 
teractions. The time dependent CP asymmetries S are shielded by large SM backgrounds 
for Bd-Taeson, whereas in Bg system the SM expectation is ~ in many cases (e.g. 
Bs J /^<P) and the unparticle scenario might reveal itself. 

We have seen that CP violation in Bd D^D^ and B ^ tv can be maximal in the 
unparticle scenario. After this phenomenological section we shall elaborate on whether 
a CP asymmetry in leptonic decays is possible. Thereafter we shall turn to the question 
of whether the scale invariance at the TeV-scale or near scale-invariance could still be 
effective at heavy flavour scales ~ 5 GeV. 



4 Constraints from CPT on (new) CP-violation 

The invariance under CPT symmetry imposes constraints on the amount of CP-violation; 
it enforces the equality of the partial sum of rates of particles and antiparticle^ to be 
made more precise below. Neither the SM nor any well-known new physics model predict 
CP-violation in leptonic decays such a.s B —>■ tu studied in this paper. The aim of this 
section is to verify explicitly whether the CP-violation is consistent with the constraints 
from CPT. 

Let us note that we expect that CPT-invariance holds for a theory with a local 



hermitian Lagrangian such as in Eq. (10). The explicit verification of CPT invariance 



demands that QC{x)Q ^ = C^{—x) = £(— x), where 9 = CPT denotes the com- 



bined CPT-transformation. The Lagrangian (10) fulfills this requirement provided that 



am grateful to Ikaros Bigi for drawing my attention to this fact. 
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QOu{x)Q^^ = Olf{—x), which we cannot verify exphcitly since we do not have equations 
of motions or a Lagrangian for the unparticle field at hand from where we would infer the 
transformation under C, P and T. There also exists a general proof of the CPT-theorem 
in the framework of axiomatic field theory [51] based on general principles and axioms 
such as Lorentz invariance, uniqueness of the vacuum and causality of field commutators. 
Concerning the latter we would like to mention that we have seen in a previous section 
that the unparticle field obeys causality, c.f. Eq. (|8]). Summarising, although we are not 
able to prove CPT-invariance we at the same time do not find any indications why it 
should be violated. 

It is well known that CPT symmetry implies equality of the decay rates of particles 
and antiparticles. In practice there is even a stronger consequence , e.g. [19], [33] or [50] 
where it was applied to charmless S-decays. The final state particles can be divided into 
subclasses of particles which rescatter into each other. It is a fact that the sum of the 
partial rates of these subclasses for a particle and its antiparticle must be the same. This 
can be inferred from the following relationship [19] between the weak decay amplitudes 
of a i?-meson and its antiparticle 5 to a final state 

{f.\H^ec.y\Br = $^(/.|5V.)(/d^dccay|5) , (45) 

i 

where -ffdecay corresponds to the weak transition operator and S is the scattering matrix. 
This relation is derived from the completeness relation 1 = \ fi){fi\ and the fact that 
the CPT-operator is antiunitary. An equivalent but alternative relation on the level of 



decay rates can be found in Ref. [33]. From Eq. (45) it is then inferred that all states fj 
which rescatter into form a subclass whose partial rates of particles and antiparticles 
sum to zero 

J2mB^f^) = 0, (/,|^t|/.)-^o (46) 



where 



AT{B^f)^T{B^f)-T{B^f) . (47) 
The exact relation between the CP asymmetry and the difference of decay rates can 



be infered from Eq. (17). Whereas the new CP asymmetry generated by Acp{D^D^) ~ 
Ar{Bd D^D^) may be compensated by AT^Bd DqDq) for instance, it is at first 
sight not clear which mode would compensate for the new CP asymmetry in Acv{t^) ~ 
AT{B^ r+z/). Among the SM final states there does not seem to be an appropriate 
candidate. We are led to look in the unparticle sector for a suitable candidate. A firm hint 
can be gained by counting the coupling constants. Denoting the weak coupling by v and 



the unparticle coupling by A (10), the CP asymmetry, which arises due to an interference 
of the two amplitudes depicted in Fig. [l| is of the order 0(A^f ^). The processes B^ —>■ 
with an interference of the two amplitudes depicted in Fig. [7] has the same counting in 
the coupling constants. One amplitude corresponds to a tree decay and the other one 
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incorporates a virtual correction due to a fermion loop of the r and the v. The process 
is kinematically allowed since the unparticle has a continuous mass spectrum. 
It does not proceed at resonance, but rather behaves like a multiparticle final state and is 
a realisation of Georgi's observation that the unparticle field in a final state behaves like 
a non-integral number du of massless particles. 





Figure 7: Decay — > lA^ ^ the double lines denote an unparticle (left) leading order (right) 
with virtual ri^-loop correction . 

We shall now explicitly verify the CPT constraint 

Ar(fi+ ^ r+z/) + Ar(fi+ ^ W+),,_ioop = . (48) 
For the sake of simplicity we shall assume as previously that there is no flavour dependent 



perturbation in the neutrino sector and that X^p = —X^s in (10). The formula for the 
flrst difference can be read off from Eq. (27) 

Ar(5+ ^ T+v) = -ABf^ sin(0) sm{du7i)Aru (49) 



• tA\ f2 

— sm(0j;^^= rrirj^ 



2\/2n rrib 



B 



m 



2 X 2 



m 



B 



-^f !/ -^ub^ub Urv\ Ad, 



■U\ A 2 



Note that the cancellation of the phase factor sm{duT^) by the same factor in the denom- 
inator, as previously mentioned, is crucial for the cancellation here since the graphs in 
Fig. [7] do not involve this factor! The amplitude of the graph in Fig. [7] to the left is 



A 



ub * 



A{B U )Fig[7l;icft) = A;^ * ^1 
and the amplitude of the graph to the right of Fig. [7] is 

AiB^ ^ W^)Fig[7];right) = ^S*^l\^TvA2 

G p 



(50) 



(51) 



where we have factored the weak parameters in A(i^2)- The fermion-loop Tls-p is given by 
the correlation function 

ns-p(p| = ml) = i ! d'xe-'^--^ (0|T[z7(l-75)r](x) [f{l-^,)v]mQ) . (52) 
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The decay rate is calculated from 



2m 



B 



d<l> 



with / = Adj^{m%Y^ 



being the phase space volume. The difference of decay rates is given by 



2 \du-2 
B) 



(53) 



(54) 



Since Ai is real only the imaginary part of A2 will enter. The only strong phase is due 
to the r and the u going on-shell in the loop in Fig. [T] (right). Therefore we only need to 
know the imaginary part of the fermion loop which is given by 



Im[ns^p(m| + zO)] 



An 



-m, 



m 



B 



(55) 



Assembling the formulae we get 



sin(0) 



2\/2'K mb 



m 



B 



I A 



'wl A 2 



/-I 



,(56) 



which fulfills the CPT constraint Eq. (48) together with (49). 



We have explicitly verfied the CPT constraint (46) for the decay B ^ rv with 
unparticle-SM interactions given by the Lagrangian (10). We do not dare to speculate in 
any detail on how a decay might be observed in a laboratory experiment. It 

can be said though that the unparticle has directed momentum, mass and charge which it 
directly inherits from the 5-meson. Moreover in the case where there is a CP asymmetry 
in 5 — > rz/ due to unparticles, it is precisely the CPT constraint (48) which tells us that 



there is an excess of charged unparticle degrees of freedom produced. Whether a part of 
this charge could annihilate into neutral particles or decay into charged particles remains 
unclear since the nature of this degree of freedom remains unknown at this stage. These 
questions could be addressed once a concrete model realising the unparticle scenario is 
known. 



5 Breaking of scale invariance - dimensional analysis 

The SM at the electroweak scale is not scale invariant. The logarithmic running and 
in particular the vacuum expectation value of the Higgs, which give masses to the fun- 
damental particles, are responsible for the breaking of scale invariance. It is therefore 
a legitimate question at what scale the symmetry breaking will be transmitted to the 



unparticle sector by the the effective Lagrangian (10). This will depend on the strength 
of the coupling and the relevance of the operators in the latter. 

The authors of reference [S] have addressed this question, which we shall adapt accord- 
ingly for the weak sector. Assuming the unparticle field couples to an operator acquiring 
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a definite mass scale, the latter will break scale invar iance at some energy. Let us assume 
for instance that the Higgs couples to the unparticle operator, following Ref . [5] , as follows 



^" = T^l^^rO„. (57) 
with Xh = c^{Au/Muy^^o~^ in our notation. We have used a new symbol Ouq for the 



unparticle operator. This operator is not the same as the one used in Eq. (10) since it 
has to be electrically neutral. The important question for the analysis in this paragraph 
is what the value the anomalous dimension du assumes. In the case where we think of 
the unparticle as being charged under SU{2)l, Ouo would appear as SC^^ ~ 'i{.l^)(lOuo 



in addition to the effective Lagrangian (10) and dy^^ = dy seems unavoidable. In the case 
where Ou is the only unparticle field then Ouq = OuO^ would be a composite field with 
anomalous dimension in the range < du^ < 2(i^, where the Thirring model at coupling 
A = 27r [55] would be an example saturating the lower bound and supersymmetric QCD at 
the conformal IR fixpoint [52] an example saturating the upper bound. In the following 
we shall quote values for the bounds and the mean value explicitly. The Higgs VEV 
= is expected to break scale invariance at a scale A 



^^yK^u,=x' K = Ku{\H^y~'-o . (58) 

What would this scale be in the cases we have investigated? Besides Ku there are two 
unknowns in the equation above, first du which appears explicitly in our results and 
Xh = {An / Mi(Y'^^o~'^ , In the latter the matching coefficient will remain unknown but 
we can extract the ratio {Ku/Mu) from pdd{Ptu) in terms of the UV dimensions. Taking 
B — »• DD as an example the breaking scale is 

where Rbd = \c's'c'^s\/\^chVcd\ = 1, deviating it from 1 corresponds to a readjustment of 
Pdd in terms of the matching coefficients. Assuming for example Au = ITeV, du = 1.2, 
duo = (0, 1.2,2.4), pdd = 10^^'^, the ratio of amplitudes and the breaking scale for fixed 
values of UV dimensions become 



Add ^ -0.40i?DD 

{duv,duvo) = (3,6) A ^ (67,20,1) GeV(i?DDC^)'/('-°''-'''-'^ 

{duv, duvo) = (3, 3) A ^ (300, 180, 50) GeV(i?-^/'c^)i/(^-°'2-«'i-6) . (60) 

The situation is not conclusive, which is not surprising bearing in mind that in the absence 
of a model there are simply to many unknowns. In the case where both UV dimensions 
are the same, which should be the case when Ou and Ouq result from the same structure. 
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a small matching coefficient Cy is needed for a sizable effect at the heavy flavour scales. If 
the UV dimensions differ by a factor of two, which is the case when Ouo = OuOjj, effects 
are possible for moderate matching coefficient . 

The effect Add = —0.40 appears larger than the analysis or conclusions in Ref. [5] 
suggest. There are two reasons. First and simply, the CP violating phenomenon investi- 
gated in this paper is linear in the ratio of amplitudes, whereas [5] describes a case where 
the effect is proportional to the square of the amplitude. Secondly it was assumed that 
the SM Lagrangian has dimension four. The crucial point is that the weak Lagrangian 
has dimension six, dc^^^y. = 6 being suppressed by two powers of the weak scale, whereas 
the unparticle Lagrangian has dimension dc^^^ = djj + dsu- In terms of the effective 
Lagrangian (10) and the Eq. ([9]), 4 < dc^^^ < 5 the unparticle operator is more relevant 



than the weak operator. This gives rise to an enhancement factor in the amplitudes 



which is explicit in the results of Eqs. (21), (41). In more physical terms one could state 



that the weak boson propagates at the high weak scale whereas the unparticle propagates 
at the low heavy flavour scale. 

Adapting the analysis of Ref. [S] we imagine an experiment at a scale /iRF, the un- 
particle Lagrangian (10) scales as Ccs = Xs / A^~^^^^^^ '^Vhf*''''^"! the weak Lagrangian as 
-^weak — GpfJ'UF and the ratio is 

Imposing that the energy scale of the experiment is higher than the breaking scale, i.e. 
HuF > A , the following bound is obtainec'^' 



This equation is easily interpreted. The first factor measures the ratio of the two couplings. 
The second is a measure between the relevance or dimension of the SM operator that is 
coupled to the unparticle and the dimension of the Higgs operator. In the third term the 
scale of the experiment has to compete with the Higgs VEV. The fourth term is peculiar 
to the weak interactions, as described above, and is due to the fact that the weak process 
takes place at the weak scale Gp^ and the unparticle propagates at the low scale fip. 
The fifth term is due to the difference of anomalous dimensions of the charged unparticle 



operator in the effective Lagrangian ( 10 ) and the neutral unparticle operator coupling to 



the Higgs VEV (57), whether it acts as an enhancing or decreasing factor depends on the 



anomalous dimensions. In a later paper [H] further dimensional analysis is explored. It 



^Setting cu ^ 1, the fourth and the fifth term to one and taking the square root of the equation, the 
bound in Ref. [Sj is recovered with = e. 
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is observed that when the couphng A multiphed by the supression scale i/A^""^*^'^^^"^-' is 
combined into a single scale then under the assumption 1 < d^^ < 2 < duv, 

it is inferred that A2 < Mu < ■ ■ ■ < A4 < A3 which seems counterintuitive at first sight 
since higher dimensional operators could receive an enhancement. 

Dimensional analysis is not very reliable. The construction of an explicit model would 
help to answer questions and presumably constrain the structure of the effective La- 
grangian (10). 

Possible candidates are extensions of the standard model featuring near conformal 
dynamics, such as the walking technicolor theories. Those theories are close to an infrared 
fixpoint and hence have slow varying coupling constants. A complete extension of the SM 
featuring walking dynamics and its link to the underlying gauge theory has been given in 
Ref. m. 



6 Critical discussion and conclusions 

In this paper we have investigated the consequences of the unparticle scenario in heavy 
fiavour physics. The new feature is a CP odd or strong phase that arises in the prop- 
agator as a consequence of the non integral scaling dimension. This gives rise to very 
characteristic and novel CP violating phenomena. 

The drawbacks of the scenario are that there is as yet no concrete model and that 
it is not clear to what energies the scale invariant sector extends. The lack of a model 



is overcome by parametrizing an effective Lagrangian, c.f. (10), at the cost of many 
unknown coefficients which have to be constrained. We have investigated the extension 
of the scale invariant sector to lower energies resorting to dimensional analysis. We have 
found that effects at the heavy flavour scales are possible provided the coupling of the 
unparticle field to the Higgs VEV is moderate at the scale Ay. The effects are sizable for 
two reasons. Firstly the scaling dimension of the unparticle Lagrangian is more relevant 
than the one of the weak Lagrangian and secondly the effect of CP violation is linear and 
not shielded by a large SM background. 

Bearing in mind the breaking of scale invariance we have chosen decays where the 
unparticle propagates at a relatively large scale. The two examples we have investigated 
are the decays —>■ r^u and — »• D^D~. In doing so we have assumed the scale 
invariant sector extends to the scale ~ 5 GeV for the former and to ~ 2 GeV for the 
latter. We have not considered decays into final state particles as in for instance Ref. [Tj 
They would also lead to signals but we have assumed the unparticles to be weakly coupled. 

We have chosen cases where the SM is described by a single weak amplitude and the 
unparticles add a second weak amplitude with strong phase allowing for the CP violation. 
In this sense our analysis does not differ from other model analyses with two amplitudes. 
The particularity of the unparticle scenario as compared to other models is that it is 
an example where the large strong phase might be generated by the strength of the 
coupling constant and that the contribution to other (fiavour)-channels is qualitatively 
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different from otlier models, for example from those generating the strong phases through 
penguins. 

The prediction of a CP asymmetry Acp = —C for leptonic decays seems a unique 
feature of the unparticle model, which has puzzling consequences to be discussed below. 
The reparametrization invariant is the product of two quadratic invariants, one from 
the quark sector and one from the lepton sector. As an example we have looked at 
T^v in conjunction with the constraints from the branching ratio. Generic plots 
for the parameter space of the anomalous dimension and the weak phase difference are 
shown in Fig. |3j Maximal CP violation is possible for certain values of the parameter 
space. The current experimental data is not yet strong enough to set absolute bounds. 
Comments on flavour related decays are given in 3^ In particular the channel D ^ nu 
might be of interest since more events have been collected [39j than m B ^ tv [37] , . 
To the knowledge of the author there are no experimental data available with bounds on 
CP asymmetries in leptonic decays. Charge symmetry is usually implied in the analysis. 

The investigation of the non- leptonic decay Bd D^D^ was motivated by the large 
asymmetry Cd+d- reported by Belle [43]. We have neglected the penguin contribution 
and treated the decay in naive factorization. As compared to B —>■ ru there is a third 
observable, the time dependent CP asymmetry Sd+d-. The latter agrees rather well with 
the SM predictions and sets constraints on Cd+d-. It is possible though to find values 
where the CP violation is maximal and satisfies the constraints of the branching ratio 
and the time dependent CP asymmetry. As for B ru, plots for generic parameters 
are shown in Fig. |6} It is encouraging that for small ratios of effective couplings the 
constraints from Sd+d- allow for a large negative asymmetry Co+d- as reported by Belle 
whereas the opposite sign seems to be disfavoured. This fact is general to any analysis 
with two amplitudes as outlined in section 2^ the unparticles just provide a scenario 
with two amplitudes and possible large weak and strong phase differences. The true 
meaning is that in the case where the decay is described by two amplitudes, the sign of 
the Belle measurement is more consistent than the opposite sign. Discussions on U-spin 



and colour related decays are given in section 3.4 Let us emphasize two points from 



this section once more. Generically we would expect a large asymmetry in Cd+d- to be 
accompanied by a large asymmetry in the color related Cj/^ttq- Currently the experimental 
value C^f^^^^ = —0.11(20) [35] is not conclusive and moreover B^ — » J/^^ttq and on the 
theoretical side, complications arise due to non-factorizable contributions. For Bd decays 
the time dependent asymmetries are typically proportional to sin(2/?) or sin(2a), the large 
angles of the Bd triangle, and new physics contributions are therefore hard to see. For Bg 
decays, the mixing phase is 0^ ~ and therefore the unparticle scenario could give rise 
to sizable corrections. This would be particularly interesting for Bg —>■ J /'^(j) which aims 
at the extraction of the Bg mixing phase 0^ at the LHCb. 

We have verified in section |4] that the novel CP-violation satisfies constraints from 
CPT-invarince, namely the equality of the sum of partial rates, of the subclasses of final 
states rescattering into each other, of particle and antiparticle. Since the SM and no 
well-known new physics model predicts a CP asymmetry for leptonic decays such as 
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T^v we have inferred that the compensating mode must be due to unparticles. 
As we have quantitatively verified, the compensating mode is — >■ lA^ . This might 
appear surprising at first sight but is possible since the unparticle does not have a definite 
mass but a continuous spectrum like a multiparticle state which was one of the basic 
observations in Georgi's first paper [1]. 

Clearly the unparticle scenario would benefit largely from the construction of an ex- 



plicit model. The question of the breaking of scale invariance and what a rea unparticle 
in a laboratory experiment would mean could be addressed and it would presumably also 
provide structural constraints on the coefficients of the effective Lagrangian. 
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A Some plots as a function of du 




'As opposed to a virtual particle, on which we focused throughout this paper. 
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Figure 9: The ratio of unparticle to SM model amplitude A^-j^pD) in units of the ratio of 
effective couplings [Pt-!^(dd) " 10^] (30)(44), versus the anomalous dimension du. (left) A^u (30), 



(right) Add (44) 



B Explcit results in coordinate space 

B.l The commuator 

The commutator of the unparticle field 

C{x) = {0\[Ou{x),OuiO)m (A.l) 
may be obtained from the time ordered product 

Cpix) = {0\TOuix)OuiOm (A.2) 

via the general formula 

C{x) = -2isgn{xo)lm[CF{x)] . (A.3) 

The correlation function Cf{x) is obtained by Fourier transformation of ([6]), c.f. Ref. |i53j 

Cf{x)= [ ^e-^-(-zA^(P^)) = - ^'7' o .^^'l ^ J}"^""] ^ (-^' + ^0)-'" • (A.4) 



(27r)4 ' '^^ 7r2 2sm{dun)T{2- du) 

The imaginary part is 

lm[{-x^ + iOy^^] = - sin(rfi^7r)e(x2)(x2)-'^" (A.5) 
and we obtain the commutator 

C{x) = -zsgn{xo)e{x'){x'r'- sinidun) ^ . ^ 



isgn{xo)Q{x'){x\''-{du - l)^?|r^^^|^^-r^ (A.6) 



-l/2-2dw 

V{2du)T{2 

The free field case du ^ 



lim C{x) = — sgn(xo)(5(x^) (A. 7) 

du^i+ 2n 



may be recovered by use of the formula lime^o+ e^l^ = 6{z). Or for any integer n 

lim C{x) ~ -isgn{xo)6^'"'~^\x^) , 

it is seen that the commutator has support on the light-cone only. 
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B.2 The Thirring model - an example with phase factor 

The Thirring model belongs to the class of exactly solvable two dimensional models. It 
is a fermionic model with a vector current-current interaction. The exact solution of the 
two point function was obtained by Johnson [54J as a function of free fields 

{0\T^{x)^{0)\0) = -ze-'^'^^^^'^Gola;) , (A.8) 

where 

Doix) = ^ log{-x' + zO) Goix) = (A.9) 

are the free bosonic and fermionic Greens functions. We have identified 7 = (A^/47r^)(l — 
A^/47r^)~^, where A is the current-current coupling constant, as for instance in [55]. N.B. 
7 > in accordance with ^ and du = I + 'J- We recover 

(0|T*W*(0)|0> = ^(3^^^ (A.10) 



the formula (A. 4) in the fermionic case up to an overall normalization. The phase factor 
arises due to resummation of thresholds at > 0. Note that the overall normalization 
in a scale invariant theory is a matter of convention and is hidden in the arbitrary scale 
factor in the logarithm of the free bosonic function log[(— -|- iO)/i^]. In the notation 
of Eq. [2] the scale /i is proportional to the fixed point scale Ay. This scales exhibits the 
phenomenon of dimensional transmutation. 

C Final state interaction in Bd D^D~ 
consistent with naive factorization 



In this appendix we shall obtain the isospin final state interaction phases within the naive 
factorization approach. The isospin analysis from K ^ tttt and B vrvr, c.f. [33] is 
transferable to Bd D^D~ [H]. The D-mesons are 1 = 1/2 states. Angular momentum 
conservation implies that only 7 = and 1 = 1 states are formed as final states in the 
decay. Denoting the amplitudes 

A+- = A{Bd D+D-) A+- = A{Bd D'^D'') = A{Bu D+D^^) , 

isospin symmetry implies 

A+' = A, A+- = ^(A + Ao) AO° = ^(A-Ao) (A.ll) 
from where the famous isospin triangle follows 

A+o = A+-+A°° . (A. 12) 
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Let us introduce the following notation 



(A.13) 



where we have factorized the final state interaction phase in the corresponding isospin 
channels. 

In the naive factorization approximation Zq = Zi = Z and therefore 

A+o = Zcosf'^^')e^(^-^«)/2 



,Zsin l^.l^)e^i^^-^o)/2 



^00 

= Ze'^' . 

The isospin triangle becomes rectangular 

Two out of the three rates have been measured [35] 

B{Bd D+D-) = 1.9(6) • 10-^ B{B+ D^D+) = 4.8(1) • 10" 



Neglecting irrelevant phase space effects we obtain from (A. 14) 

- 0.4(2) . 



cos 



5i - 5c 
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